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Abstract 

We present a method to compute, by numerical simulations of lattice QCD, the 
light-cone wave functions which enter exclusive processes at large momentum transfer, 
such as electromagnetic elastic scatterings, or exclusive semideptonic decays as B —* n 
(B — ► p) and radiative decays as B — > K*j. The method is based on first principles 
and does not require any model assumption. 



In this paper we propose a method to compute, by numerical simulations on the 
lattice, the light-cone wave functions which allow to predict the form factors relevant 
in many exclusive processes, such as electromagnetic elastic scattering at large mo- 
mentum transfer, or exclusive semi-leptonic decays as B —> tt (B — > p) and B — > if* 7 
decays [jl]]-|§- Our approach allows to compute the form factors in exclusive semi- 
leptonic B decays at small values of the invariant mass of the lepton pair, q 2 , a region 
which is not accessible by standard lattice techniques, which are confined in the region 



- 2 ^ 

'lmax 



~~ <r hi -| 11]. We also show that, from suitable combinations of lattice correla- 



tion functions, it is possible to determine the full light-cone wave function, denoted 



generically as <&, and not only its moments, as done in the past [12]. The possibility 
of computing directly $ is a significative advantage: higher moments of the light-cone 
wave functions, being related to higher dimension operators, are in general afflicted by 
power divergencies when a hard cutoff [] is used (or renormalon ambiguities in dimen- 
sional regularizations) . This makes very problematic the definition of the renormalized 
operators, i.e. those which have finite matrix elements when the cutoff is removed fl3fl . 
With our method, instead, no renormalization is needed as we get directly <3? from the 
appropriate correlation functions. The technique described below strictly follows a sim- 
ilar proposal recently made to determine the shape function for inclusive heavy-hadron 
decays 



We now explain how the method works. The light-cone wave functions are universal 
quantities, i.e. quantities independent from the process at hand (electromagnetic elastic 
scattering at large momentum transfer, exclusive B decays, etc.). Thus, in order to 
illustrate our proposal, we start by considering a very simple prototype which enters, 
for example, the elastic j*(q) + vr — > tt scattering Q and B — ► tt semileptonic decays || 

q) = if d 4 x e^ x (7r-(p n )\d(x)^ l5 S(x; 0)u(0)|0) (1) 

where, for simplicity, and without loss of generality, we have taken S(x; 0) to be the 
scalar Feynman propagator, satisfying the differential equation 

- D 2 S(x;0) = 5 4 {x -0) . (2) 

With this choice the multilocal operator appearing in eq. (|l]) is gauge- invariant. 
F^(Pn, q) can be written as 

F^p^q) =ijd*x (7T-(p;)|d» 7 ^ 75 5Q(x;0)n(0)|0) , (3) 
In the lattice regularization the hard cutoff is given by the inverse lattice spacing I /a. 
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where Sq(x;0) = e iq ' x S(x;0). For simplicity we consider first the case of a pion with 
a small momentum (at rest), i.e. \p n \ <C \q\ {p-w = 0). For — q 2 = Q 2 3> p 2 = M 2 , 
we can separate the large frequency modes (~ Q = \J Q 2 ) from the low energy modes 
(~ Aqcd) and expand Sq(x;0) in powers of Aq C d/Q 

0) - ( n2 . ,* p>m- ) ' ^ 
\-Q z + 2iq- D + ieJ (x . 0) 

In eq. (|j) we kept only the leading terms of the expansion in powers of 1/Q, and those 
which become leading when u = Q 2 /{2p 1T ■ q) ~ 1. Using eq. (Q) one finds 

q) = M C du n ®: {u) — (5) 

Jo -Q z + 2uq ■ p w + le 
where the light-cone wave function $> n (u) is defined through the relations 

(7r-(p n )\d(0h^ 5 (iD^)...(iD^)u(0)\U) = -iUM n p^ ...p^ 

+ iB n 5^p^ + ... , (6) 

with 

M n = [ duu n $> n {u). (7) 
Jo 

In eq. @, the contributions proportional to B n are suppressed by higher powers of 
1/Q in all the relevant kinematical region. 

We now consider the Fourier transform of F^(p n ,q) defined as 

F"(t,p„q) = J ^e-^F^q) . (8) 
For t > 0, by closing the contour of integration over go below the real axis, we find 

F»{t,p«,q) = -iUp% / du = = ■ . (9) 

Jo q 2 + u l E 2 + 2uq ■ p n 

If, for consistency with the order at which we are working, we neglect the terms of 
0(M 2 /Q 2 ) and 0(pi/f), we get 

rl -j(-M_B 7r + -y/ q 1 +2uq-p 7r )t 

F»(t,&,$) ~ -i/.K / du —====- . (10) 

•'0 2^q 2 + 2uq-p 7r 

It is convenient to write eq. (^) as follows 

F^t,p n ,q) = -if v p£ / du „ F ^ , , (11) 

Jo 2uE 7T (p 7T + q/u) 
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where E n (p) = \Jp 2 + M 2 . In the latter form, one recognizes that eq. (||) is valid for 
arbitrary pion momenta. In the general case, eqs. (p|)-(|TTD can be derived from the 
following expression 

q) = UP% t du w2-u->(*i U \ ' (12) 

Jo W 2 + 2(q + p T ) ■ k + le 

where W 2 = (q + p n ) 2 = Q 2 (l — x)/x is the squared invariant mass of the recoiling 
hadron system and k = up w — p^ is the residual momentum of the struck parton 
((D^ 1 . . . D^ n ) ~ k^ 1 . . . k^ n = 0(Aq CD )). The expression above shows that different 
dynamics which occur depending on the value of x. At small x, we can retain only the 
W 2 term in the denominator and the amplitude becomes proportional to the lowest 
moment of fJv which is unity; when W 2 ~ Q^-qcd, ah the moments become comparable 
and we need the full wave function finally in the region where W 2 ~ Aq CD , the 
light-cone approach fails. This is the analog of what happens in inclusive decays at the 



end point of the energy spectrum [14]. 



It is straightforward to make the analytic continuation of the above expression to 
the Euclidean space-time which is used in numerical simulations 

F»(t,p n ,q) = -iUp£F(t,&,$ 

= -iUft / du , (13) 

Jo 2uE n (p 7T + q/u) 

By studying the time-dependence of F^it^p^, q) at several values of q and p n , we can 
unfold the integral above and and extract the light-cone wave function <fr w . 

We give below the practical recipe to implement the calculation of F^lt,^,^) in 
lattice simulations. From the definition of F^ip^, q) in eq. (|l]), one has 

F»(t,p w ,q) =ijd 3 x e-*'* (Tr(p w )\d(x,t) 7 ^ 5 S(x,t-,d,0)u(0)\0) . (14) 

The above amplitude can be extracted directly from a suitable ratio of lattice three- 
and two-point correlation functions 

Rfr&rf = Hm e- E ^%^ P :f , (15) 



where E w = \/M 2 + p 2 ; 

G%(t f ,t,p n ,q) = I d 3 x e^ 3 (0\U^(t f )d(x,t)^ 75 S(x,t;6,0)u(0)\0) (10) 



and 

= (0111^(^)^(0)10) . (17) 
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Hp v (t) and At$(t) are the pion interpolating field and the axial current (A tl (x) = 
u{x)^^^d{x)) with definite spatial momentum p n 

D^(i) = J <Pxj*'-%A' l (2,t) , A^(t) = J d 3 x^ s A"(x,t) . (18) 

In the Euclidean, using the transfer matrix formalism, we have 

n^(f,t) = e 6t U^(x)e- 6t , (19) 

so that the correlation functions have an exponential dependence on the energy of the 
external states. This implies that, in the limit tf — > oo, the lightest state, corresponding 
to a pion, dominates the correlation functions ( |16| ) and (|T7|), since all higher-energy 
states are exponentially suppressed. In this limit, with t > 0, we then obtain 

R(t,p w ,q)^ F(t,p n ,q), (20) 

which is the desired quantity. 

It is obvious that the same technique can be used to compute all the light-cone 
functions which have been introduced in the literature. In the case of the pion, for 
example, those denoted as & p and 3> CT in ref. |^] correspond to the non-local amplitude 
of eq. (Q), where 7^75 is replaced by 275 and a^ u respectively. The same holds true, 
taking into account all the complications entailed by the polarization of the physical 
states, in the case of the 3>s which are relevant for the p [||], the K* meson 0] and the 
nucleon || cases. 

In practical calculations one may use either the scalar propagator defined in eq. @) 
or the approximate propagator introduced in eq. (|j). In the general case, the latter 
can be written as 

= ( ? I WT^T-Z \ n^- ) (21) 
\{q + Pn) 2 + MQ + Pn) -D + ieJ 



(x;0) 



Sq(x) can be written as 



e +iQx+/2 



Sq(x) = S LEET {x) , (22) 

where x + = n ■ x, with n M = (q^ + p%)/Q, and S LEET (x) is the light-cone propagator 
of the Large Energy Effective Theory (LEET) [jlEj, which satisfies the equation 

2iD + S LEET {x) = 5\x) . (23) 

is the appropriate vector which becomes light-like in the elastic region, n 2 = 
W 2 /Q 2 ~ Kqcd/Q when W ~ Q^-qcd- Thus, the extraction of the <I>s from Sq(x) is 
equivalent to the use of the LEET. Note that the calculation of the physical light-cone 
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function ^ w using the full propagator of eq. (|2|) does not need any renormalization 
(for <E>p and <3? CT , the same renormalization constants of the operators tpj^ip and ifja^ip 
must be applied). The same calculation using the S LEET requires, instead an overall 
(further) logarithmic renormalization of the amplitude [|16| , which can be computed in 
lattice perturbation theory. The ultraviolet divergences of the LEET correspond per- 



turbatively to infrared divergences in the full theory [17]. In the latter case the infrared 
divergences are automatically regularized by the non-perturbative contributions in the 
physical matrix elements and no renormalization is required. 

As discussed in ref. |l4||, it is not clear whether the use of S LEET (x) will be convenient 
in practice, since this propagator is much more singular than the full one. For this 
reason we expect that the correlation functions computed in numerical simulations 
using S LEET (x) will be affected by larger statistical fluctuations. 
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